Graphene edge structures: Folding, scrolling, tubing, rippling and twisting 
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Conventional three-dimensional crystal lattices are terminated by surfaces, which can demonstrate 
complex rebonding and rehybridisation, localised strain and dislocation formation. Two dimensional 
crystal lattices, of which graphene is the archetype, are terminated by lines. The additional available 
dimension at such interfaces opens up a range of new topological interface possibilities. We show 
that graphene sheet edges can adopt a range of topological distortions depending on their nature. 
Rehybridisation, local bond reordering, chemical functionalisation with bulky, charged, or multi- 
functional groups can lead to edge buckling to relieve strain, folding, rolling and even tube formation. 
We discuss the topological possibilities at a 2D graphene edge, and under what circumstances we 
expect different edge topologies to occur. Density functional calculations are used to explore in 
more depth different graphene edge types. 



INTRODUCTION 



II. TOPOLOGICAL DISTORTIONS WITHIN 
THE GRAPHENE PLANE: FLAT EDGES 



A finite material is necessarily terminated by an inter- 
face. While the bulk of a crystalline material respects 
crystal symmetry, this symmetry is broken at the in- 
terface. Material interfaces are thus heterogenous and 
generally more reactive than the material bulk. Symme- 
try breaking can lead to imbalanced local strain which 
needs compensating in some way, typically through in- 
terface relaxation but also potentially through disloca- 
tion creation. Interfaces also create dangling bonds and 
enhanced reactivity, which once again can be mitigated 
through various effects such as chemical rehybridisation, 
interface reconstruction, and chemical functionalisation. 

In the case of three-dimensional bulk solids their 
interfaces are two dimensional surfaces, but for two- 
dimensional materials such as graphene, their interfaces 
are one-dimensional lines. As we discuss in this article, 
this difference in dimensionality means that graphene in- 
terfaces have a number of possible topological distortions 
and relaxation modes which are not available at 'classic' 
surfaces of three-dimensional materials. The result is a 
rich variety of potential interface types in graphene, all of 
which can radically alter the properties of the material, 
even at quite long range from the interface. 

In the following article we discuss different interface 
behaviour in graphene. We start by discussing more clas- 
sic interface behaviours, giving examples of rehybridisa- 
tion, interface reconstruction, restructuring and chemical 
functionalisation. We then consider key new topologi- 
cal distortions that exploit the additional available third 
dimension. For this it is useful to consider a geometri- 
cally anisotropic sample of graphene such as a graphene 
nanoribbon. We take here the example of a graphene 
ribbon of infinite length, i.e. with two principle axes, 
one along the ribbon length and one orthogonal to it, 
both of which can demonstrate distinct topological edge 
distortions. 



By fiat edges we imply restricting the edge to remain 
within the graphene plane. In this case edge effects are 
direct one-dimensional analogues of surface behaviour in 
three-dimensional crystals. 



A. Rehybridisation 

Cutting carbon along the armchair, or 'boat', direc- 
tion, [1010] results in a series of atom pairs. These 
are able to rehybridise from sp 2 in the graphene bulk, 
towards sp triple bonding to stabilise the edge. This 
shortens the bond from 1.412 A in the graphene 
bulk to 1.22 A, creating one of the most stable non- 
functionalised graphene edges, with a formation energy 
of only 1.10 eV/Ai, as shown in Figure [TJleft). 



B. Interface reconstruction 

Dangling bonds on a two dimensional surface can be 
saturated through local rebonding between atoms at the 
interface. This results in a surface superlattice containing 
more than one bulk lattice vector in at least one direction; 
for example the (10 0) surface of Si consists of Si atoms 
with three neighbours which stabilise by bonding in pairs, 
creating a [2 x 1] surface reconstruction. 

A direct analogy at a graphene edge can be observed 
for the so-called Klein edge 2 , cutting graphene along the 
[2110] direction in the shuffle plane. The result is a line 
of chemically unstable singly coordinated carbon atoms 
with correspondingly high edge formation energy (2.21 
eV/Ai). Allowing symmetry breaking along the edge 
direction results in pairwise reconstruction, giving the 
graphene equivalent of a [2 x 1] surface reconstruction (a 
[2] edge reconstruction) which drops the edge formation 
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Figure 1: Different edge structures for flat unterminated graphene edge (adapted from Supplementary Materials, Reference 
Ivanovskaya et al. 1 ). (left) Armchair edge showing rehybridisation, (center) Zig-zag edge showing restructuring, with alternate 
edge bonds rotated about their bond centers, and (right) Klein edge showing edge reconstruction through pairwise rebonding. 
Edge formation energies (eV/A) are indicated below each structure. 



energy to 1.50 eV/A 1 , as shown in Figure Upright). 



(a) 



C. Interface restructuring 

In some ways a subset of interface reconstruction, 
this class of edges undergo extensive reordering in the 
interface layer in order to create a more energetically 
favourable surface structure. Cutting graphene once 
again along the [2110] direction, but this time in the 
glide plane, results in a zig-zag terminated edge, slightly 
more stable than the reconstructed Klein discussed above 
(1.34eV/Ai). This edge is metallic and cannot stabilise 
through local rehybridisation. However a series of 90° 
bond rotations along the edge around the bond centres 
for alternate bonds changes the crystal lattice from a 
hexagonal array to a periodic array of alternating pen- 
tagons and hexagons 3 , a linear analogue of the proposed 
layered Haeckelite structures. In this case the edge atoms 
are now arranged pairwise, similarly to the armchair 
edges discussed above, and the edge carbon atoms can 
thus once again rehybridise, resulting in a very stable 
edge structure with formation energy of only 1.09eV/A\ 
as shown in Figure Reenter). 



III. TOPOLOGICAL DISTORTIONS 
ORTHOGONAL TO THE RIBBON AXIS: 
FOLDING, SCROLLING AND TUBING 

We now consider possible edge distortions which have 
no direct analogy in surface termination of three dimen- 
sional materials, i.e. these interface modes are unique to 
two dimensional layered materials. 

If we consider the graphene nanoribbon in cross-section 




Figure 2: Schematic cross-section showing four possible out- 
of-plane distortions of graphene sheet edges orthogonal to the 
edge line: (a) flat, (b) folded, (c) folded and rebonded (tubing) 
and (d) scrolled edges. 



looking along the axis, there are four distinct classes of 
edge distortion possible. The first is to simply remain 
planar, as discussed above (Figure [2^). However if we 
consider a thought experiment where the edge is now 
pulled up out of the graphene plane and back above the 
graphene layer, there are three possible structures de- 
pending on the angle at which the edge then approaches 
the graphene surface below it. If the approach angle is 
less than 90° the edge folds back on itself, resulting in 
a bilayer structure (Figure ^jp). If the approach angle is 
90° the edge will bond into the sheet below (Figured), 
creating a pseudo-nanotube at the graphene edge, and 
if the approach angle is greater than 90° the edge folds 
back in on itself, creating a scroll (Figure [2jl) . We now 
discuss each of these in more detail. 
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A. Folding 

The flexibility of graphene layers promotes the for- 
mation of self- folded nanostructures (Figure [2)3). High 
resolution transmission electron microscopy (HRTEM) 
studies 4 — have shown that graphene edges can fold back 
on themselves. The energetic cost of bending the layer 
is compensated by Van der Waals interactions in the 
stacked region 12 . Multiple folding of graphene results 
in multi-layer regions and highly curved folding edges 13 . 
The critical self-folded length at which point folding be- 
comes thermodynamically favourable has been theoreti- 
cally deduced 14 . 

Graphene folding can occur in any direction, however 
the interaction between the folded layers depends on 
the resultant atomic stacking. Scanning tunneling mi- 
croscopy (STM) and HRTEM studies suggest the possi- 
bility of both AA and AB stacking for folded graphene 
layers^ 5 -. Energetic comparison for folded structures 
demonstrates the preference of graphene folding: the 
global minimum is associated with AB stacking of the 
entire flat region, while the local minimum with a mix- 
ture of AB and incommensurate stacking occurs in the 
presence of a small twist of folded graphen o 16 i 17 . 

Folded graphene edges present a combination of a 
nanotube-like and multilayer graphene structures which 
give rise to peculiar electronic propertie d 13 ! 16 . Nonethe- 
less such edges are potentially very stable and commonly 
seen in multi-layer graphitic stacked materials. In such 
cases it is also common to observe more than one sheet 
folded simultaneously together resulting in a bi- or multi- 
layer fold. 

B. Scrolling 

Folding and consequential rolling up of a graphene 
layer into a spiral structure at the sheet edge leads to 
the formation of a scroll (Figure [2]i). A large variety 
of possible scroll structures can be obtained by coiling a 
single or multiple graphene sheets, changing the number 
of coils and sliding relatively adjacent layers. The open 
and highly modulable structure of scrolls suggests poten- 
tial applications for hydrogen storage or for use as ion 
channels^ Experimentally, scrolls have been obtained 
via arc discharge, chemical treatment of graphite or 
graphene but an easy and reproducible route for scroll 
synthesis still remains to be developed. 

Several theoretical works have investigated the forma- 
tion and structural stability of graphene scrol b 19 i 2Q . Sim- 
ilarly to graphene folding, the formation of scrolls is dom- 
inated by two major energy contributions: the energy 
increase due to the bending of the planar layer and the 
energy gain due to the van der Waals interactions be- 
tween the rolled layers. Beyond a critical diameter value 
these scrolled structures can be energetically more stable 
than the equivalent planar configurations 2 ^^, however 
in order to obtain a scroll a large energy barrier due to 



the bending rigidity should be overcome. 

Scroll formation occurs spontaneously when a critical 
overlap between layers is achieved for the partially curled 
sheet, and interestingly scroll unwinding has been ob- 
served during charge injection^ 2 -. Thus, electrostatic con- 
trol of the wrapping appears feasible, opening the way to 
possible technological applications 23 . 

The minimum innermost radius of nanoscrolls was also 
calculate d 22 ! 24 ! 25 , and is experimentally observed to be 
in the range 20-50 A 24 . Thus scrolling will be limited at 
graphene edges of larger flakes and will not be important 
for graphene nanoribbons. 

Similarly to nanotubes, the electronic structure of 
nanoscrolls has been shown to depend on their chiral- 
ity (n,m) 24 . Armchair nanoscrolls were found to be 
metallic or semimetallic depending on their sizes. Metal- 
lic scrolls have a larger density of states at the Fermi 
level than metallic single- walled nanotubes. Zigzag nano- 
scrolls were found to be semiconducting with energy gap 
much smaller than corresponding zigzag nanotubes. The 
optical properties of carbon nanoscrolls have been stud- 
ied as well: the calculated reflection spectra and loss 
function showed features of both single- wall carbon and 
multiwall carbon nanotubes 25 . 



C. Tubing 

By folding an unterminated graphene edge back on it- 
self and rebonding it into the graphene layer it is possi- 
ble to create 'nanotube-terminated' edges 1 (Figure [2^). 
Depending on the edge chirality, zigzag and armchair 
edges lead to the formation of armchair and zigzag tube- 
terminated edges, respectively. The edge dangling bonds 
are thus replaced with sp 3 -hybridised carbon atoms. For 
sufficiently large tubes (armchair tubes larger than (8,8) 
and zigzag larger than (14,0)) this results in lower for- 
mation energies than any other non-functionalised edge 
structure discussed here. 

Simulated high resolution trasmission electron mi- 
croscopy images are found to be similar to those of free 
standing edges, the primary difference being minor vari- 
ations in image contrast, suggesting 'tubed' edges may 
be difficult to distinguish experimentally. 

The electronic properties of such tube terminated 
edges present an interesting combination of graphene and 
nanotube behaviour. Rolled zigzag edges serve as metal- 
lic conduction channels, separated from the neighbouring 
bulk graphene by a chain of insulating sp 3 -carbon atoms, 
and introduce Van Hove singularities into the graphene 
density of states (Figure [3]). They may provide a way to 
stabilise and protect from chemical attack the disperse 
Fermi level state seen along metallic zigzag edges. 1 
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Figure 3: 'Tubed' graphene edge, where a non-functionalised 
edge is rolled back on itself and bonds into the layer below, 
creating a line of sp 3 -hybridised carbon atoms (marked in 
yellow), in this case for a (8,0) zig-zag 'tube', (right) Density 
of states for a (8,8) armchair 'tube' edge, where the strong 
Fermi level peak corresponding to metallic states along the 
carbon atoms neighbouring the sp 3 -bonded atoms is clearly 
visible. Figure adapted from Reference Ivanovskaya et al.— . 



IV. TOPOLOGICAL DISTORTIONS PARALLEL 
TO THE RIBBON AXIS: RIPPLING AND 
TWISTING 

A. Rippling via functionalisation 

Pristine graphene edges have dangling bonds at the 
edge atoms. Simple H-termination is the simplest way 
to saturate these dangling bonds 26 . With this approach 
very little strain is induced to the edge. Adding differ- 
ent atoms to the edge such as -F, -CI or more complex 
functional groups such as -OH or -SH changes this simple 
picture. In general most functional groups add significant 
strain along the ribbon edge, through steric hindrance, 
electrostatic repulsion between groups, inter-group bond- 
ing (such as hydrogen bonding), etc. 

This strain is energetically unfavourable, and can be 
relieved via out-of-plane distortion^. Specifically, hy- 
droxyl (-OH) terminated graphene nanoribbons of dif- 
ferent widths (notably armchair edges) have been shown 
to compensate the induced strain by forming a localised 
static ripple along the edge2£. This rippled edge is more 
stable than any fiat configuration. The strain is relieved 
via distorting the edge carbon hexagons pairwise up and 
down periodically. The ripple is localised at the edge, 
confined by the sp 2 tt - system of graphene trying to stay 
flat. These ripples are of a different length scale and 
underlying mechanism to previously observed statistical 
thermal fluctuations, and bending of the graphene sur- 
face via the physisorption of molecules such as H2O— »2S. 

In Figure [3^ relaxed armchair graphene nanoribbon 
of width 7 with thiol functional groups is shown (width 
nomenclature is taken from Reference 30 ). A key con- 
sequence of these functionalised nanoribbon edges is 
that both electronic and mechanical properties can be 
tuned. By changing the functional groups the band gap 
can vary for certain ribbon widths up to 50 %. Rippled 
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Figure 4: -SH terminated armchair graphene nanoribbon of 
width 7. (a) Rippled graphene ribbon edges with grafted SH- 
groups (in gray the graphene ribbon network, in blue the -SH 
groups, one -SH group per carbon edge atom), (b) as guide 
for the eye the effective rippled armchair graphene nanoribbon 
edges are shown by a thick red line. 

structures can also decrease the Young's Modulus for 
small width graphene nanoribbons up to 40 % compared 
to H-terminated armchair graphene nanoribbons 2 ^. 

In general the influence of "complex" functionalised 
graphene edges become increasingly important as the sur- 
face area to edge length ratio decreases 2 ^. Thus edge 
rippling is likely to be of key importance in graphene 
nanoribbons with small widths (< 2 nm) and small di- 
ameter graphene flakes. 

B. Twisting via functionalisation 

An alternative way to relieve edge strain proposed re- 
cently is twisting of the whole graphene nanoribbon—. A 
schematic is shown in Figure [5j In this study a -F ter- 
minated armchair graphene nanoribbon of width 7 was 
twisted helically and found to be most stable with a twist- 
ing angle of 4.2 per unit cell. In this case the sp 2 tt - 
system of the graphene nanoribbon is slightly distorted 
throughout compared to a fiat graphene sheet. 

The literature study used tight binding and small cells 
with twisted boundary conditions, and it is not clear the 
cell size was sufficient 31 . We have therefore also calcu- 
lated the structure and energy of such twisted ribbons us- 
ing a large orthorhombic unit cell containing 1620 atoms 
(ribbon width 7) performing a 360° twist in the unit cell 
with an angle of 4° between fundamental unit cells, us- 
ing a more accurate DFT-LDA method as implemented 
in the AIMPRO code^ 2 -). Such system sizes are possible 
due to development of a new filtration method for DFT 
calculations (for computational details se o 33 i 34 ). First 
preliminary results suggest that twisted ribbons are less 
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Figure 5: Schematic armchair graphene nanoribbon width 7, 
twisted via 30 ° between each original unit cell. The red and 
green line are guides for the eye indicating the two ribbon 
edges. This twist is exaggerated, previous calculations find 
twist angles closer to 4° (see text) 

stable than the flat ribbons with rippled edges. This 
seems sensible since edge rippling disturbs only the edge 
of the sp 2 -graphene network whereas twisting distorts the 
entire ribbon. Even if twisting can occur this suggests a 
limitation in width where twisting can be applied, as rip- 
pling the edge is independent of the graphene nanoribbon 
width. 



V. CONCLUSION 

We have discussed here a range of topological possi- 
bilities for graphene edges, many of which are unique 
to two-dimensional crystal lattices. In-plane possibilities 
include rehybridisation, restucturing and reconstruction, 
all of which have been observed in HRTEM. Out-of-plane 



distortions depend on the axis of distortion and include 
folding, scrolling, 'tubing', as well as rippling and possi- 
bly even twisting once chemical functionalisation of the 
graphene edge sites is included. Indeed it should also 
be possible to combine these orthogonal distortions, re- 
sulting in, for example, rippled edges which then form 
scrolls. 

We note that the discussion here applies to free- 
standing graphene. Many topological distortions rely on 
a balance between an energetic cost to distort the sheet, 
offset by either strain relief or inter-layer interaction en- 
ergy. Adding a further surface energy via interaction 
with a substrate will alter these ratios, and at the very 
least change the zones of stability for each topological 
distortion. The study of interfaces in two-dimensional 
materials opens up a rich diversity of structures which 
we have only begun to characterise and exploit experi- 
mentally. 
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